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We report on a 59 Co NMR investigation of the trigonal cobaltate CasCo2 06 carried out on a single 
crystal, providing precise determinations of the electric field gradient and chemical shift tensors, and 
of the internal magnetic fields at the non-magnetic Col sites, unavailable from former studies on 
powders. The magnetic field-induced ferri- and ferromagnetic phases at intermediate temperature 
(e.g. 10 K) are identified by distinct internal fields, well accounted for by purely dipolar interactions. 
The vanishing transferred hyperfine field at the Co I site indicates that the Co (I) orbitals do not 
participate in the intra-chain superexchange, in disagreement with a previous theoretical model. 

The strong Ising character of the system is confirmed experimentally by the field dependence 
of the resonance lines, indicating that local moments are saturated even at the phase boundaries. 
In the vicinity of the critical fields, nuclear spin-spin relaxations detect the spin reversal dynamics 
of large magnetic assemblies, either Ising chain fragments or finite size domains, which drive the 
metamagnetic transitions. Such collective excitations exhibit a glassy behavior, slowing down to 
subacoustic frequencies and freezing at low temperature. The relevance of such slow fluctuation 
modes for the low-temperature multi-step behavior reported in the magnetization is discussed. 

PACS numbers: 75.30.Kz, 75.30.Et, 75.60.Jk, 76.60.-k. 



I. INTRODUCTION 

The nature of magnetic order in the geometrically frus- 
trated trigonal cobaltate Ca3Co20g, an Ising-type mag- 
net exhibiting low-dimensional behavior, has been a puz- 
zle for over a decade. The crystal structure of Ca3Co206 
(Fig. [1} may be viewed as the result of the alternate 
stacking along the c axis of face-sharing CoOg octahedra 
(Col) and Co06 trigonal prisms (Coll) with a rather 
short Col - Coll distance c/4 = 2.59 A, building up 
chains arranged over the ab plane in a triangular lat- 
tice, with a much larger spacing (inter-chain distance 
a/V3 = 5.24 A) ±£ The Co 3+ oxidation state at both 
sites is now established^"^ contrary to earlier predictions 
for a Co 2+ -Co 4+ charge disproportionation£ 

Magnetism resides on the Coll site, where a trigo- 
nal crystal field yields a high-spin ground state with 
a large unquenched orbital component ((L z ) w 1.7, 
(2S Z + L z ) = 5.3), while Col is in a low-spin state 
(i.e. non-magnetic) £ Spin-orbit coupling gives rise to 
a single- ion anisotropy term — Z?5 2 of order 70 meV 
(corresponding to an anisotropy field B anis sa 200 T 
along c), whence the strong Ising-type character of mag- 
netism in this system^" 6 - Exchange coupling between 
the Co 3+ II ions is ferromagnetic (FM) and relatively 
strong along the chains, whereas it is antiferromagnetic 
(AF) and much weaker between chains. While the dom- 
inant intra-chain exchange interaction promotes quasi 
one-dimensional Ising-type ferromagnetism, the actual 
magnetic structure is the result of the competing residual 
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FIG. 1: (Color online) Lattice structure of Ca3Co206 in the 
hexagonal setting. Left: perspective view of the spin chains 
formed by the CoOs polyhedra along the c-axis; right: cross- 
section in the ab plane, evidencing the arrangement of chains 
in the triangular lattice. The dark and light polyhedra repre- 
sent the CoOe octahedra and C0O6 trigonal prisms, respec- 
tively. 



inter-chain AF coupling, giving rise to a complex mag- 
netic phenomenology as a function of temperature and 
magnetic field. 

Susceptibility measurements actually demonstrate the 
onset of magnetic order at T < T c « 25 K.— ^ The nature 
of magnetic ordering however, has been the subject of 
a long-lasting controversy. Denoting spin-up, spin-down, 
or magnetically disordered Ising chains as f , 1 , 0, respec- 
tively, theoretical models indicated either a ferrimagnetic 
(FI) [ttl] or a partially disordered antiferromagnetic 
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(PDA) structure [ j" I 0] as the stable zero- field arrange- 
ments of the three inequivalent chains on the triangu- 
lar lattice in the presence of geometrically frustrated AF 
inter-chain coupling^— Early neutron diffraction experi- 
ments seemingly supported either a PDA— or a Fl£ or- 
der, although poor quantitative agreement of both struc- 
tures with magnetic peak intensities was later pointed 
outpii Recent x-r&yH and neutron scattering^ 3 - data on 
single crystals have revealed the incommensurate charac- 
ter of the magnetic order in this compound. They show 
that the actual structure, referred to as modulated PDA 
(MPDA), rather consists of long- wavelength longitudinal 
spin-density modulations along each chain, phase-shifted 
by ±27r/3 with respect to the neighboring chains Nev- 
ertheless, the MPDA phase is apparently metastable be- 
low 14 K, as indicated by the decrease in the intensity of 
the MPDA magnetic peaks.— ~— 

The nature of the magnetic order in an externally 
applied magnetic field is also far from clear. Increas- 
ing magnetic fields applied along c tend to stabilize 
the FI phase and eventually FM inter-chain order. At 
5K < T < T c , the initial saturation of the magneti- 
zation to a magnetic moment Mpi sa 1.7 hb per for- 
mula unit at /iqH > B c \ ss 0.1 T, followed by a step to 
Mfm ~ 3Mpi at ijlqH > B C 2~3.5T, actually points to 
perfect FI and FM phases below and above H C 2, respec- 
tively. However, the M(H) curves at lower temperatures 
display a puzzling multi-step behavior reminiscent of the 
quantum tunneling of magnetization detected in molecu- 
lar magnetsi 14 ! 15 A number of complex magnetic super- 
structures, distinct from both the FI and FM spin-chain 
arrangements, have been proposed in order to explain the 
fractional magnetization values. However, such super- 
structures cannot account for all the plateaus observed 
in M(H)M- 

In this paper, we address the issue of the field- 
dependent magnetic order of Ca3Co206 by means of 59 Co 
nuclear magnetic resonance (NMR). Previous 59 Co NMR 
research on Ca3Co206 polycrystalline samples, though 
successful in demonstrating the presence of non-magnetic 
Co 3+ , could only provide very limited information on the 
magnetic structure, due to the angular averaging of the 
complex multi-line quadrupolar spectra into broad pow- 
der patterns^ The FM and FI surroundings, on the con- 
trary, can be clearly resolved in our NMR investigation of 
a high-quality single crystal, while the two field-induced 
transitions, namely MPDA-FI and FI-FM, were precisely 
identified at 10 K through step-like jumps in the internal 
field and peculiar features of the spin-spin relaxations, 
indicative of glassy spin dynamics. Our study also yields 
a detailed characterization of the chemical shift and the 
electric field gradient (EFG) tensors at the non-magnetic 
Col site. 

The paper is organized as follows. Sample preparation 
and the experimental methods are briefly described in 
section [TTJ NMR spectra, providing information on the 
magnetic structures, and nuclear relaxation data, prob- 
ing the dynamical excitations of the spin system, are pre- 



sented in sections IlII Al and IIII Bl respectively. The cal- 
culated dipolar fields at the Col sites in the FI and FM 
phases are compared with experimental values in section 
HVl The results are discussed in section [V] Two appen- 
dices review in some detail the complex dependence of 
the 59 Co resonances on magnetic and quadrupolar inter- 
actions, and the longitudinal relaxation function of a 7/2 
nuclear spin. 



II. EXPERIMENT 

The investigated sample was a single crystal of size 
1.2 x 1.2 x 5.8 mm 3 (the long side parallel to the c axis) 
taken from a batch of several needle-shaped crystals ob- 
tained by the flux method. The crystals were grown using 
the following procedure: a mixture of Ca3Co20g powder 
and KCO3 (the flux), in a weight ratio 1/7, was heated 
up to 990° C in an alumina crucible for one hour and then 
slowly cooled down to room temperature. A sample of 
the same batch was emp loyed for the resonant x-ray scat- 
tering study of Ref . [l2|. The high quality of the crystals 
was confirmed by preliminary x-ray diffraction, energy 
dispersive x-ray, magnetization, and specific heat mea- 
surements. 

Magnetization measurements were performed after 
zero field cooling, by using a vibrating sample magne- 
tometer (Oxford Instruments) equipped with a 12 T mag- 
net, with a field sweep rate of 1 T/min in order to mini- 
mize any relaxation effect. Susceptibility measurements 
were carried out by means of a superconducting quantum 
interference device magnetometer (SQUID) with mag- 
netic fields up to 5T and temperatures down to 1.8 K 
(Quantum Design MPMS). The magnetic measurements 
were carried out for the two geometries corresponding to 
the magnetic field applied parallel or perpendicular to 
the c axis, i.e. parallel or perpendicular the direction of 
the chains and the spins. In the case of parallel geom- 
etry the crystals were aligned using the magnetic field 
itself at 5 T as explained in Ref. [l7|- The alignment of 
the crystals in the perpendicular geometry can be easily 
obtained thanks to their rod-like shape. 

The NMR experiments were performed by means of a 
home-built phase-coherent pulsed spectrometer— and a 
fast-sweeping cold-bore cryomagnet (Oxford Instruments 
EXA) equipped with a variable temperature insert as a 
sample environment. The crystal was mounted on a sam- 
ple rotator, providing an angular span of ±120° with re- 
spect to the applied field. Misalignment of the rotation 
axis (nominally perpendicular to the field) was estimated 
not to exceed 5 degrees. In the case of the c axis parallel 
to the field, however, a precise alignment of the sample 
was achieved, thanks to its large magnetic anisotropy, by 
loosely mounting it in the coil. 

NMR spectra were detected by exciting spin echoes by 
means of a standard P — r — P sequence, with equal rf 
pulses P of duration 3-5 us and intensity suitably reg- 
ulated for optimum signal, and delays r kept as short 
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as possible with respect to the dead time of the res- 
onant probchead («10-25/zs depending on the working 
frequency). Recording was carried out either by tuning 
the probehead at discrete frequencies in a constant field 
(frequency sweep mode), or by varying the applied field 
and exciting the resonance at a constant frequency (field 
sweep mode). Although an indirect method, we preferen- 
tially employed the latter whenever possible (namely, far 
from the metamagnetic transitions) as it yields smoother 
data, independent of the frequency response of the spec- 
trometer, in a fully automated procedure. 

Nuclear spin-spin relaxations were determined by the 
decay of the signal amplitude in the same spin echo se- 
quence P — t — P, as a function of variable delay r. Spin- 
lattice relaxations were measured by the signal recovery 
following a so-called fast saturation of the observed nu- 
clear transition (see Appendix IB]) . obtained with an ape- 
riodic train of 10-15 pulses. 



III. EXPERIMENTAL RESULTS 




A. NMR spectra 
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FIG. 2: (Color online) Frequency-swept 59 Co spectrum at 
10 K from the majority spin-up chains of the FI phase, 
recorded in a field of 0.38 T applied along the c axis. Solid 
lines are bimodal fits of the satellite peaks, vertical dash- 
dotted lines mark their centers of gravity. 



We report separately the Co resonance spectra 
recorded in the longitudinal geometry (external field 
B ext || c axis), a geometry that allowed us study the field- 
induced magnetic order in the system, and in transverse 
applied fields. Interpretation of the spectra is relatively 
simple and model-independent in the longitudinal case, 
due to the collinearity of all the fields (internal and ex- 
ternal magnetic fields, and the EFG). In the general case, 
the dependence of the 59 Co resonances on magnetic and 
quadrupolar interactions is rather complex, as recalled in 
some detail in Appendix [AJ 

1. Longitudinal applied fields 

A 59 Co NMR signal could be detected below 15 K in 
moderate longitudinal external magnetic fields B ext > 
0.1 T. Such field values assign the system to the FI phase. 
For such field values the FI phase is induced in the sys- 
tem. In contrast, no signal could be detected in the 
MPDA phase at B ext < 0.1 T. 

A typical frequency-sweep spectrum, recorded at 10 K 
in 0.38 T, is plotted in FigJ5J The spectrum consists 
in a septet of quadrupole-split Zeeman transitions v n , 
n = —3, ...,n = 3 (Eqn. I A2[) . with constant frequency 
spacing f„+i — v n = 2.11(1) MHz, in agreement with 
Ref. [l9|. In this geometry (9 = 0), and in the presence 
of a cylindrical EFG (77 = 0), as is actually the case 
(see Sec. IIII A 2)) . such a line spacing coincides with the 
quadrupolar coupling parameter vq (Eqn. [A3]) . The reso- 
nance frequency i^q of the central line, which is unaffected 
by the quadrupole interaction, corresponds to a field 
B nuc — 2irv / 59 j at the nucleus which is larger than B ext 
by approximately 1.2 T (here, 59 7 = 2ttx 10.10 MHz/T22 
is the gyromagnetic ratio of 59 Co). Such a field offset, 



though sizable, is nevertheless much smaller than the typ- 
ical hyperfine fields detected in magnetic Co 3+ , which are 
of the order of -10 T/^s and +60 T/jUb for the spin and 
orbital hyperfine components, respectively! 21 ' 22 This con- 
firms that the present 59 Co signals are from Co 3+ ions 
in a low-spin state, in agreement with previous NMR 
reports.— ^ We note that the dipolar field at site I gen- 
erated by two nearest-neighbor Coll moments of 5.2 [Ib 
aligned parallel to B ext accounts for both the positive 
sign of the field offset B nuc — B ex t, and its magnitude to 
within a few percent accuracy. This is the first indication 
that the internal field Bi nt ~ B nuc — B ext arises from the 
surrounding magnetic ions essentially through the dipo- 
lar coupling, as actually proven in detail by the internal 
field calculations of Sec. lIVI Given its dominant on-chain 
dipolar origin, the positive sign of Bi nt then indicates 
that the nearest-neighbor Co 3+ spins are aligned along 
the direction of the applied field. The present positive- 
offset line multiplet therefore originates from nuclei at the 
Col sites of spin chains aligned parallel to the external 
field, namely, the majority chains of the FI structure. 

Qualitatively similar positive-offset spectra could be 
reproduced in increasing external fields parallel to the c 
axis (8 = 0) at 10 K, up to B ext = 3.5 T. The resonance 
frequency of the central line follows a linear dependence 
on B ext (Fig.©, 

v = aB ext + b (1) 

with best-fit parameters a t = 10.21(1) MHz/T, b t = 
12.02(2) MHz (hereafter, f I hi the subscripts denote 
quantities related to spin-up majority and spin-down mi- 
nority chains of the FI phase, respectively). On the other 



4 



80 



70 



■n 60 



o 

> 50 



40 



30 




4 5 
B «|P1 



FIG. 3: (Color online) Central resonance frequency vo of the 
positive-offset spectra from the spin-up chains in either the FI 
(Bext < 3.5 T) or the FM phase (B ext > 3.5 T), as a function 
of B ext || c. The dashed lines are fits to straight lines for the 
two phases. The inset is a blow-up of the transition region. 
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lB nuc = 59 7 (1 + K c )(B ext + B mt ) (2) 



as detailed in Appendix [5] If the staggered magneti- 
zation is saturated all over the field span, the internal 
field Bi n t is independent of B ext and, by symmetry, also 
parallel to c. Working under this hypothesis, whose va- 
lidity is demonstrated in the discussion, a comparison of 
Eqs. Q] and [2] simply relates the coefficients a, b to Bi n t 
and the c axis component of the chemical shift tensor K c : 
a = 59 7(1 + A" C ), b = 59 j{l + K c )B mt , respectively. Best- 
fit values for a-)- and b^- then yield Bi nt ^ = 1.178(2) T, 
and a sizable K c « 0.01. 

Resonances from the minority spin-down chains are 
marked by a negative field offset, comparable in abso- 
lute value to the internal field in the majority chains. 
The signal could only be found in longitudinal fields ex- 
ceeding 2 T, since for smaller external fields its resonance 
frequency Vq rj 59 7 |i?int| ~ |-B«Bt|| falls below 10 MHz, 
outside the frequency range of our apparatus. A field 
sweep spectrum at 10 K, recorded at fixed frequency 
of 16 MHz, is shown in Fig. [4j It consists of a mul- 
tiplet of quadrupole-spit lines centered at -B„,— similar 
to those detected for the majority chains, with a cen- 
tral resonance at Bq = 2.792 T. The two satellites at 
£?_2 = 3.172 T and B-% ps 3.36 T are however reduced in 
amplitude by approximately 75% and 90% with respect 
to the symmetric peaks at B2, B3, respectively The sig- 
nal loss at B ext > B-2 is due to the approach of the 
FI-FM transition at B C 2 ps 3.5 T, which is accompanied 
by a notable relaxation phenomenon that partially wipes 
out the signal intensity over a wide field interval close 
to B C 2 (see below). The rightmost satellite at B-% is 



40 
35 
_30 
■Q 25 
-g 20 

< 

10 
5 




v L =16MHz 



2 
1.5 
1 

0.5 




1 

obi 


c 


rt 





\ « 












°c°° 


s 


a* 


* 





3.35 



3.4 



L. jL jL Jk. Jt _L _jk 



2.2 



2.4 



2.6 2.8 

B e Xt n 



3.2 



3.4 



FIG. 4: (Color online) Field-swept spectrum (B ext || c) at 10 K 
from the minority spin-down chains, recorded at a constant 
frequency v — 16 MHz. The solid line is a bimodal fit to a 
septet of quadrupolar satellites, with positions constrained by 
Eqn. IA2I The inset is a blow-up of the rightmost satellite at 
B- S w 3.36 T. 



furthermore displaced from the expected position by ap- 
prox. -2mT (figure inset), a value well accounted for by 
the demagnetization field, since the macroscopic mag- 
netization significantly deviates from the FI plateau at 
H = /i " 1 i3_3 (Fig. [TTj) . The quadrupole coupling pa- 
rameter is estimated as Vq = 1.97(1) MHz from the line 
spacing, vq = 59 7(1 + K c )(B n -i - B n ), and is slightly, 
but significantly smaller than in the spin-up chains. 

The resonance frequency of the central line follows a 
similar linear dependence on B ext (Fig. [5]), with slope 
a_i = 10.22(1) MHz/T, and a negative intercept b± = 
-12.48(3) MHz (Eqn. [J). The former coincides within 
error with a-j- , whence the same chemical shift value K c sa 
0.01 is obtained. From bj,, the internal field at Col in the 
minority chains is estimated as Bi nt i = —1.221(4) T. 

The negative-offset resonance lines from the minority 
spin-down chains completely disappear in external fields 
Bext > 3.5 T, while only a positive-offset septet can be 
detected, that is qualitatively similar to the one from 
the majority FI chains below 3.5 T and with an identical 
quadrupole splitting v Q = 2.11(1) MHz (Fig. EJ). The 
integrated amplitudes of the positive-offset spectra above 
and below 3.5 T, and proportional to the number of 59 Co 
nuclei excited in spin- up chains, have a ratio very close 
to 3:2. These facts indicate that all the spin chains are 
now aligned parallel to the external field, i.e. the system 
is in a FM state. 

The positive-offset septet from the spin-up chains could 
be followed at 10 K while increasing B ext in steps of ap- 
prox. 0.1 T across the field- induced FI-FM transition, 
without losing the signal. The resonance frequency of 
the central line is plotted vs. B ext in Fig. [31 Expcri- 
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FIG. 5: (Color online) Central resonance frequency of the 
spectra from the minority spin-down FI chains, as a function 
of Bext || c. The solid line is a fit to Eqn.[T] 



mental points clearly follow two different straight lines 
above and below 3.5 T (Eqn.[T|), which signals the occur- 
rence of a transition at a critical field B c2 = 3.50(5) T. 
Line parameters on the FM side are fitted as a FAI = 
10.23(1) MHz/T, in good agreement with the a values in 
the FI phase, and a reduced b FM = 11.40(1) MHz, cor- 
responding to a step-like decrease of the internal field, 
B™ = 1.1 15(2) T. 

The transition is also accompanied by an abrupt broad- 
ening by a factor of about 2 of the satellite lines above 
B C 2- A closer inspection of the resonance lines reveals, 
on both sides of the transition, an asymmetrical bimodal 
shape with a sharper peak and a broader negatively- 
shifted shoulder, that is more pronounced in the FM 
phase. Line shapes and widths, in contrast, are inde- 
pendent of satellite order m. These facts indicate that 
the observed inhomogeneous broadening is magnetic in 
origin. The increase in line width IS.vt, above B C 2 quali- 
tatively agrees with a dominant line broadening from the 
demagnetization field, which is proportional to its net 
magnetization, and inhomogeneous throughout the sam- 
ple due to its non-ellipsoidal shape. The agreement is not 
quantitative, however, since a step in Ai>l by a factor of 
3 across B c2 is expected from the corresponding step in 
the macroscopic magnetization. Such a discrepancy sug- 
gests the presence of an extra broadening mechanism, 
relatively more important in the FI phase. 

2. Transverse applied fields 

The crystal was also mounted with its c axis perpen- 
dicular to the external field, in order to study the trans- 
verse components of the EFG and K tensors in an az- 
imuthal scan of B ext in the ab plane. A typical field- 
sweep spectrum at 10 K, recorded at a fixed frequency 
vi, = 75.6 MHz, is shown in Fig. [7] The spectrum con- 
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FIG. 6: (Color online) Field-swept (B ext \\ c) 59 Co spectrum 
in the FM phase at 10 K, recorded at a constant frequency 
v = 64 MHz. Solid lines are bimodal fits of the satellite peaks, 
vertical dash-dotted lines mark their centers of gravity. 



sists of two septets of relatively sharp quadrupole-split 
transitions, with positive (Bq < vl/ 5 7) and negative 
offsets (Bo > vl/ 59 j) of their respective central lines. 
The integrated amplitudes of the positive- and negative- 
offset multiplets are in a ratio very close to 2. Following 
the above arguments, the two multiplets originate from 
the majority spin-up and minority spin-down chains, re- 
spectively, in the presence of a sample misalignment of 
the order of a few degrees, yielding a longitudinal field 
component of order a few hundreds mT, strong enough 
to induce FI magnetic order in the sampled The effect 
of sample rotation around an axis nominally coincident 
with the crystallographic c axis is summarized in Fig. El 
showing the dependence on the azimuthal angle </> of the 
central lines -Bnfi Bqi, and of the first-order quadrupole 
splittings, calculated as (B-\ — B\)/22& Overlaid on the 
experimental points, the figure also shows for comparison 
simulated shifts (Fig. Ek) and splittings (Fig. Eb) pre- 
dicted by Eqs. [AD and [AD for 6 = 81° (the angle between 
c and B nuc sa Bi oc , resulting from a vector addition of 
Bext = 7.45 T with a perpendicular B int « 1.2 T), and 
rhombicity factors 77 of order 0.1. Clearly, the experimen- 
tal splittings do not follow the predicted angular periodic- 
ity of 7r, rather showing a 27r period (Fig. EK)- Moreover, 
the quadrupolar shift is negligible even for large 77 values 
(Fig.Eb)j since it is a second order term in vq/vz- There- 
fore, the entire azimuthal dependence of both quantities 
must be due to a slight misalignment between c and the 
rotation axis, as well as between the latter and the plane 
perpendicular to B ext , indicating that 77 is vanishingly 
small. 

The spectrum of Fig. [7] is fitted to two septets of 
Lorentzian lines, with constrained positions calculated as 
functions of external and internal fields, chemical shift, 
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FIG. 7: (Color online) Field-swept spectrum at 10 K, recorded 
at a constant frequency of 75.6 MHz, in a nominally trans- 
verse geometry {B ex t -L c). The solid line is a global fit to 
two septets of Lorentzian lines with positions constrained by 
Eqn. IA2I Vertical dashed and dash-dotted lines mark peak 
positions for the majority and minority spin chain signal, re- 
spectively. 



quadrupole interaction, and a geometric angle a between 
the c axis and B ext . In the fit, the in-plane chemical shift 
tensor component K a , a, and vq\, vq^ (which are in prin- 
ciple independent for the two chain types) were treated 
as free variational parameters, while Bi nt ^, Bi nt i, and K c 
were kept fixed to the values determined by NMR with 
c || B ext (longitudinal geometry). A further constant 
parameter, p = dtp/dB ext = 3.4(2) x 1CP 3 T _1 , account- 
ing for a field-dependent tilting angle ip of the electronic 
moments, was calculated from the a&-plane macroscopic 
susceptibility, determined by SQUID magnctometry as 
X± = 1.8(1) x 10"Vs/T. The fit yields a negative 
K a = -6.7(2) xlCT 3 , and VQr = v Qi = 2.11(1) MHz, 
coincident with the quadrupole line splitting of the ma- 
jority signal in the longitudinal geometry, and so is in 
disagreement with the value determined for the minority 
chains in large longitudinal fields. 

B. Nuclear relaxations 

Nuclear relaxations were systematically studied on the 
central quadrupole line (m = 0) in the longitudinal ge- 
ometry (B ex t || c) at 10 K as a function of the external 
field. The longitudinal nuclear polarization recovers equi- 
librium in good agreement with Eqn. IB21 indicating that 
spin-lattice relaxations are actually driven by magnetic 
fluctuations, as expected in this system. 

Figure [9|3 displays spin-lattice and spin-spin relaxation 
rates of the positive-offset septet at 10 K over the 2-6 T 
field span, an interval comprising the FI-FM transition 
B C 2 . Spin-spin rates T 2 _1 arc larger than spin-lattice rates 
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FIG. 8: (Color online) Symbols: experimental first order 
quadrupole splitting (a) and central peak position (b), in field 
units, as a function of the azimuthal angle <j> in the transverse 
geometry. Dashed lines are guides to the eye. Solid lines: 
simulations with 9 — 81°, corresponding to B ex t J- c (see text 
and Eqn. IA2]) . for increasing 77 values. 



Tf 1 = 2W (Eqn. |B2J) by at least 3 orders of magnitude, 
and the two relaxations clearly follow quite different be- 
haviors vs. B ext . The former exhibits a peak at B C 2, 
with such a large maximum value that it nearly leads 
to the disappearance (wipcout) of the signal at this tem- 
perature. In contrast, T± exhibits a monotonic decrease 
with increasing fields, except for a shallow feature at B&. 
Moreover, the time dependence of the spin-spin relax- 
ations exhibits a marked non-exponential decay over a 
wide field interval across the FI-FM transition. The sig- 
nal decay is best fitted to a stretched exponential func- 
tion, A(t) = A exp[— (2t/T 2 )^], with stretching expo- 
nent /3 peaked to a value of 1.5 at B&, i.e. the intrinsic 
lineshape (related to the spin-spin relaxation function by 
Fourier transform) is nearly Gaussian close to the tran- 
sition. 

The difference in magnitude and field dependence of 
the Tj -1 and T 2 _1 rates is not peculiar to the longitudinal 
orientation, as it can be qualitatively reproduced with the 
c axis rotated by a — 45° in B ex t = 4-7 T (not shown), a 
field interval encompassing the FI-FM transition in that 
geometry. This indicates that the behavior of Tj~ and 
T 2 _1 vs. applied field is essentially invariant with respect 
to the orientation of the nuclear field B nuc sa _B/ oc rela- 
tive to the c axis (note that B ext 3> Bi nt in this case, so 
that the Bi oc direction is approximately that of B ext ). 

Perfectly consistent data are recorded from nuclei in 
the minority spin-down chains, for B ext approaching B C 2 
from below. A qualitatively similar behavior of the two 
kinds of relaxations is also observed at 10 K close to 
the other metamagnetic transition, i.e. from the FI to 
the MPDA phase on decreasing the longitudinally ap- 
plied B ext down to B c i sa 0.1 T (Fig. (SJi). Here again 
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FIG. 9: (Color online) Spin-lattice (triangles) and spin-spin 
relaxations (bullets) measured on the central resonance of the 
positive-offset septet at 10 K as a function of B ext || c, (a) on 
approaching the FI-MPDA transition from above, and (b) 
across the FI-FM transition. Inset of panel b: stretching 
exponent ft of the spin-spin relaxation function vs. B ex t- In 
the figure T 2 _1 rates are scaled by a multiplication factor of 
10~ 3 for clarity. 
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FIG. 10: (Color online) Spin-spin (triangles) and spin-lattice 
relaxation rates (bullets) measured in the FM phase on the 
central resonance (yo = 60 MHz, B ex t = 4.749 T parallel to 
c) as a function of temperature. The dashed line is a fit of 
T7/ 1 to an Arrhenius law. In the figure T 2 1 rates are scaled 
by a multiplication factor of 10 -4 for clarity. 



T^ 1 is orders of magnitude larger than T^ 1 and increases 
as B ext approaches the critical field, while T-j -1 is non- 
divergent. As in the vicinity of B C 2, the homogeneous 
lincshape is quasi-Gaussian, as spin-spin relaxations are 
fitted to stretched exponential decay forms, with j3 in the 
range 1.2- 1.4. 

Spin-lattice relaxations were also measured as a func- 
tion of temperature and applied field, both in the FI 
and FM phase. Results of a typical temperature scan, 
performed in a longitudinal field of 4.749 T (i.e. in the 
FM phase) are shown in Fig. [T0J The figure also shows 
the temperature dependence of T.^ 1 for comparison. The 
spin-lattice rate clearly follows a thermally activated be- 
havior T^ 1 = T cx . 1 exp— A/T, as also found by other 
authors.— The Arrhenius law is obeyed over a tempera- 
ture range of 8— 13 K at this particular field, correspond- 
ing to a T\ variation by more than two decades, with 
activation energy estimated as A = 96(2) K. Above this 
temperature interval the signal is lost due to an exceed- 
ingly short T 2 , while experimental T^ 1 points exhibit ex- 
cess relaxation at T < 8 K, indicating that the relaxation 
mechanism responsible for the Arrhenius behavior be- 
comes unimportant at low temperature, and it is shunted 
by other more effective relaxation channels. Spin-spin 
rates seemingly saturate at low temperature and also ex- 
hibit a steep temperature dependence at T > 8 K as well, 
although the evolution is clearly not Arrhenius-likc. 

Thermally activated T{ l (T) were detected at all fields, 
irrespective of whether the system exhibits FI or FM 
order. The field dependence of the activation energy 
A(B ext ), as well as its interpretation in terms of the el- 



ementary excitations of an Ising spin chain system, will 
be reported elsewhere. 



IV. DIPOLAR FIELD CALCULATION 

A calculation of the total dipolar field at the Co I site 
is required for a quantitative analysis of our experimen- 
tal internal field values. We report such a calculation 
here for each of the three magnetic environments, major- 
ity FI-, minority FI-, and FM-ordered chains. The con- 
tribution of nearby moments (referred to as the proper 
dipolar field Bdi P ) is calculated by summing up the in- 
dividual dipolar fields from Coll spins SSj at posi- 
tions r j inside a Lorcntz sphere with a radius of sev- 
eral unit cell radius centered at a Col site r$ = 0, 
B d ip = Y., ///;;;(••'>/• ,;• , • SSj - r?SSj)/r)j, while spins 
outside the Lorcntz sphere are treated as a macroscopic 
magnetic moment density, giving rise to a demagnetiza- 
tion field Bdem inside the sphere, as is customary. The 
latter is strictly constant throughout the sample volume 
only in the case of an ellipsoidal sample surface, whereby 
it is simply proportional to the saturation magnetization 
M s and to the difference in demagnetizing factors of the 
inner Lorentz sphere 1/3 and the outer sample surface 
N, Bdem = 4?r(l/3 — N)M S . For a rotational ellipsoid 
magnetized along its principal axis c, the demagnetizing 
factor N c is calculated as2£ 

Nc = MHVe) 
e = y/1 - (r a /r c y 
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TABLE I: Calculated dipolar sum Bdip, demagnetization field 
Bdem, and total dipolar field B tot = Bdi P + Bdem as a function 
of magnetic chain ordering, compared with experimental Bi nt 
values. 



Order 


Bdi P (T) 


Bdem (T) 


B tot (T) 


exp. (T) 


FM 


0.997 


0.137 


1.134 


1.115(2) 


Fit 


1.122 


0.046 


1.167 


1.178(2) 


FLj. 


-1.246 


0.046 


-1.201 


-1.221(4) 



A non-ellipsoidal sample may be approximated by an 
effective ellipsoid accounting for the mean Bd em , while 
deviations from the ideal shape gives rise to a distribu- 
tion in demagnetization fields, contributing only to the 
linewidth. We assume effective ellipsoid radii r a = r^, 
r c equal to half of the sides of our square-base paral- 
lclcpipedal specimen as a best ellipsoidal approximation. 
The choice of r a and r c does not critically affect the de- 
termination of Bdem, since A^ c depends weakly on e in 
the case of an elongated sample (e.g. dN c /de « 0.3 for 
our crystal). For instance, a 10% error on e would result 
in an error ABd em of order 5 and 2 mT in the FM and 
FI phases, respectively. 

The results of our calculations are summarized in ta- 
ble Hand compared with experimental results. Both the 
dipolar sums and the demagnetization field were cal- 
culated by assuming the Co II moments saturated to 
±g(S z )/j 1 B along c in all the three types of magnetic 
chains, with a ± sign as appropriate for the chain or- 
der, and identical absolute moment values estimated as 
g {S z } = 5.25(3) (in units of /j,b) from an extrapolation of 
M (H) data in high field. Lattice parameters were taken 
as a = 9.061 A, c = 10.367 A as reported in the litera- 
ture at 40K»i Numerical convergence of Bdi P could be 
attained with a Lorentz sphere radius of the order 30 of 
unit cells. Calculated and experimental values differ by 
less than 2% in all cases. Agreement should be considered 
satisfactory in view of the uncertainties in the magnetic 
moment (of order 1%) and the demagnetization factor. 
This confirms that the internal field is dipolar in origin, 
while transferred hypcrfine contributions, if they occur, 
are negligible. 



V. DISCUSSION AND CONCLUSIONS 

We organize the discussion of our experimental results 
into separate sections for sake of clarity, each one focusing 
on one main conclusion of the paper. 

A. Interactions of 59 Co: evidence against a direct 
Co I-Co II superexchange path 

Our NMR study on a single-crystal provides a deter- 
mination of the chemical shift and the EFG tensors at 
the Col site. The chemical shift is less precise due to 
the uncertainty in the zero-shift 59 Co reference^ and to 



the presence of a large internal field which complicates 
the unraveling of K components (see Appendix [A} . The 
absolute magnitude of K, of order 1%, though sizable 
in absolute terms, is however only slightly larger than 
the typical values reported for diamagnetic complexes of 
cobalt.— The relatively small value of K agrees with the 
known large crystal field splitting at the octahedral I site, 
of order 0.65 eV^ since a near-lying crystal field excited 
multiplct usually gives rise to a large chemical shift via 
the van-Vleck mechanism, as often encountered in non- 
magnetic transition metal ions. 

The EFG exhibits cylindrical symmetry around c, in 
accordance with the fact that c axis is a C3 symmetry 
element of the point group of Col. Consistent values of 
the quadrupolar frequency vq (proportional to the EFG) 
are estimated in different chain types and experimental 
geometries, with the notable exception of the minority 
chains in external longitudinal fields approaching the FI- 
FM transition value B&, where an EFG deviation of or- 
der -7% is found, well beyond the experimental uncer- 
tainty. We argue that a magnetoclastic coupling of the 
Co 3+ ions in counter-oriented minority chains might be 
the source of the observed EFG reduction, as already 
proposed for other cobalt-based compounds studied by 
Mossbauer spectroscopy^ 

The internal field Bi nt at the non-magnetic Col site 
can be accounted for by purely dipolar interactions with 
the surrounding Coll spins, within an absolute accuracy 
of 2 mT in all of the three chain types. This result, along 
with the experimental value of K, sets a very stringent 
upper limit to a transferred hypcrfine field at the nucleus. 

The contact hyperfine field may originate from both a 
tiny on-site magnetic moment on Col, as that invoked 
in a neutron scattering refinement,— and from the cou- 
pling to nearest neighbors Coll spins. These two con- 
tributions might partially cancel out, since either sign is 
possible.— Additional inter-chain contributions are negli- 
gible in view of the large chain spacing. Nevertheless, we 
can safely rule out an on-site term, since a residual on- 
site moment would require an admixture of the ti g with 
the excited e g orbitals, in contrast with the large crystal 
field excitation energy implied by the observed moder- 
ate chemical shift value. For instance, a huge chemical 
shift K sa 5 was reported for the resonance of 141 Pr in 
low-spin Pr 3+ , in combination with a small fractional 
moment « 0.01 hb on the Pr ion^ Therefore, only the 
transferred contribution is to be considered. 

An intra-chain transferred hypcrfine term B^f would 
offset the total dipolar field B to t (Tab. [TJ) by ±1^/1 for 
the Flf and FM chains, and by =p|B^/| for the FLJ. chains, 
depending on the relative sign of B to t and Bhf- However, 
a finite Bh / biasing term does not improve the agreement 
between the calculated B to t±Bhf and the experimentally 
determined internal fields. We conclude therefore that 
the slight discrepancy between the experimental and the 
calculated internal fields must be due to other sources 
of error, and that the contact field, if present, is much 
smaller than 2 mT. 



A negligible hypcrfinc field at the Col site is an un- 
expected result. Transferred hypcrfinc fields at nuclei 
of non-magnetic ions in magnetic compounds of tran- 
sition metals are typically of the order of several hun- 
dreds mT or more, unless perfect cancellation occurs at 
a symmetric site of an AF structure, which is not the 
case in Ca3Co206- For instance, the transferred hyper- 
fine field at 139 La in the FM manganite Lai-^Ca^MnOs 
equals 3T, i.e. approximately 10% of the on-site hyper- 
fine field at the Mn sites, where the magnetism resides^ 
Here, a virtually vanishing contact field at Col indicates 
that there is no significant hybridization between the 
Co 3+ (I) and the Co 3+ (II) wave functions either directly 
or through an oxygen bridge. 

The mechanisms underlying exchange and hypcrfinc 
transfer are similar, although not exactly identica h 31 i 32 
Therefore, the lack of a detectable transferred hypcrfinc 
field strongly suggests that Co 3+ (I) ions do not partici- 
pate in the exchange interaction J\ between two neigh- 
boring on-chain Co 3+ (II) spins. We stress that this find- 
ing constitutes experimental evidence against the cur- 
rently accepted theoretical model by Fresard et al, who 
instead proposed an exchange path across Co I assuming 
a large overlap of Co 3+ (II) with Co 3+ (I) wave functions, 
yielding a direct integral t = 1.5 eV.— 

B. Nature of the metamagnetic transitions: perfect 
Ising behavior. 

At 10 K two types of field-dependent magnetic struc- 
tures are detected, namely perfect FI and FM ones, in 
accordance with the two plateaus in magnetization data 
at the same temperature. The two structures are clearly 
identified by distinct local fields at the Co I sites in good 
quantitative agreement with the dipolar fields calculated 
for the corresponding spin arrangements. The transition 
from a FI to a FM phase shows up as a sharp step in Bi nt 
at a threshold field B c i ~ 3.5 T, revealed by an abrupt 
crossover of the resonance frequency vs. field between two 
distinct v^(B ext ) functions (Fig. [3]). 

Quantitative comparison between NMR and magneti- 
zation data in longitudinal fields at 10 K (Fig. ITTj) pro- 
vides an insight into the magnetization process. Close in- 
spection of the magnetization curves as a function of the 
applied field H = B ext j\iQ shows that the steps appear 
smoother in the M(H) curves than in the NMR frequency 
data vo(B ext ) (figure inset). Moreover, the plateaus in 
M(H) are only approximately flat, with a moderate field 
dependence that is larger in the FI state (for instance, 
H^dM/dH w 0.05/j b /T for IT < fi H < 3T). 

In principle, incomplete saturation of M(H) over the 
plateaus might be due to either unsaturated c axis com- 
ponents of individual moments in a homogeneously mag- 
netized sample (e.g. due to spin canting), or inhomo- 
geneous magnetization, although the former could not 
easily be reconciled with the known strong single ion 
anisotropy of the system. NMR demonstrates experi- 
mentally that the cobalt spins are collinear and saturated 




FIG. 11: (Color online) Magnetization curves at 10 K as a 
function of the external field H applied along c. Inset: detail 
of M(H) across the FI-FM transition, overlaid for compari- 
son to the frequency offsets Auq = i/q — aB ext of the central 
NMR line of either Fit or FM chains, where a is the best-fit 
parameter to Eqn. [1] (whence Afo ~ b tx Bi nt within error). 

even in the vicinity of the FI-FM transition. Proof is pro- 
vided by the field proportionality coefficient a of Eqn. [TJ 
whose experimental values for the three chain types were 
found to coincide within an error Sa w 10 kHz/T. The ob- 
served value of dAI/dH contributes a 15 kHz/T slope (i.e. 
of order Sa) to v>o(B ext ) via the demagnetization field, 
regardless of the mechanism responsible for it, and inde- 
pendent of the magnetic order in a chain. A hypothetical 
field dependent canting, however, would affect the dipo- 
lar field at Col, which is larger than the demagnetization 
field by an order of magnitude (Tab. UJ. Unless we as- 
sume equal canting in the FM, FI spin-up, and the FI 
spin-down chains, which is clearly unphysical, the lat- 
ter would produce differences in the a coefficients of the 
three lines well beyond experimental errors, and contrary 
to experimental evidence. 

As seen from NMR, the system apparently exhibits 
perfect Ising behavior. Therefore, the rounded steps in 
the M(H) data at the field- induced metamagnetic tran- 
sitions, with unsaturated magnetization over broad H 
intervals near them, must be due to a fraction of misori- 
ented domains or spin chains. 

C. Two distinct magnetic excitations probed by 

Tf 1 and T^ 1 

Further insight on the spin reversal process close to 
the MPDA-FI and FI-FM transitions is provided by the 
large and divergent spin-spin relaxations, accompanied 
by much smaller and non-divergent spin-lattice relax- 
ations. We recall that the spin-lattice relaxation rate 
T^ 1 probes transverse random magnetic fields fluctuat- 
ing at the Larmor frequency lol , while the spin-spin rate 
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T-2 is the sum of a term proportional to (the popu- 
lation term) and a secular term arising from longitudinal 
random fields fluctuating at virtually zero frequency . 34 ! 35 
In the presence of isotropic and fast fluctuations in the 
so-called narrowing limit (such that uil is much smaller 
than the roll-off frequency of the fluctuation spectrum) 
the two rates tend to coincided Here, the similar behav- 
ior of relaxations measured in longitudinal and oblique 
external fields rules out the possibility that anisotropic 
fluctuations, namely, entirely along the c axis and hence 
ineffective for spin-lattice relaxations, are the source of 
the huge difference in T\ and X2 time scales. Excess spin- 
spin relaxation and T 2 _1 peaks must instead depend on 
fluctuations in a regime opposite to the narrowing limit, 
namely with characteristic frequencies much lower than 
ujl and hence contributing only to the secular term in 
T 2 -1 . Such a regime, which has already been reported in 
the literature for other oxides,— is referred to hereafter 
as the nearly-static limit. 

To this end, the Kubo-Tomita theory of line shape 
appropriate for continuous- wave magnetic resonance 3 ^ is 
modified in the context of pulsed NMR by an experimen- 
tal low-frequency cutoff, namely the reciprocal time win- 
dow uj c = (2t)~ 1 of a spin-echo experiment (here, r is the 
time separation of the two excitation pulses). Random 
fields fluctuating at frequencies u <C w c ~ 2x 10 4 s _1 may 
be viewed as static over the excitation-detection transient 
and contribute only to the inhomogeneous line broaden- 
ing, which is refocused by the spin echo sequence. Con- 
versely, fluctuations at w c < u < effectively relax the 
spin-echo signal. 

Within the nearly-static limit, the decrease of T 2 _1 with 
decreasing temperature (Fig. [T0|) then indicates the col- 
lapse of the fluctuation spectrum below the cutoff uj c at 
lower temperatures, down to a complete freezing of the 
spin dynamics observed at T < 5 K. It is worth noting 
that evidence for comparable time scales in the spin dy- 
namics was also reported from x' an d x" ac magnetic 
susceptibilities in the zero-field phase, both showing a 
strong frequency dependence over the acoustic and suba- 
coustic ranged The nearly Gaussian T2 relaxation form 
is another signature of nuclear spin depolarization from 
quasi static random fields^ In our case, random fields 
are clearly electronic in origin, since the calculated sec- 
ond moment of the nuclear dipolar fields yields a T 2 _1 
rate smaller than the observed ones by at least one order 
of magnitude. 

The peculiar behavior of T^ 1 rates therefore reflects 
the exceedingly slow spin dynamics, as in the case of 
massive excitation modes involving the collective motion 
of several spins. In view of the strong Ising character 
of the system, such excitations probably consist in the 
coherent reversal of large portions of the spin chains or 
domains. The slowly fluctuating random fields respon- 
sible for spin-spin relaxation are then straightforwardly 
identified with the stray dipolar fields originating from 
such transient defects in the magnetic structure. Also the 
peaks in T 2 ~ 1 and the stretched exponent j3 at the meta- 



magnetic transitions finds a natural explanation within 
this framework. 

A peak in T 2 _1 vs. B ext may arise in principle from 
either a field-dependent correlation time of fluctuations, 
or a peak in the random field amplitude. In the nearly- 
static fluctuation regime, however, field dependence of 
the correlation time would imply faster spin fluctuations 
at the peak. The latter is apparently inconsistent with 
the peaks in /? vs. B ext reaching nearly- Gaussian values 
of 1.5 at B c i and B C 2, a signature of a nearly frozen spin 
dynamics, following the above argument. The T 2 _1 peak 
must therefore originate from a corresponding increase of 
the quasi-static magnetic disorder, namely a larger num- 
ber of misoriented domains or chain fragments, close to 
the transitions. We note that the peaks of f3 vs. B ext 
coherently reflect changes in the spatial distribution of 
random fields on approaching the metamagnetic transi- 
tions. According to a well-known mechanism , 34 i 40 the 
transition from a quasi- Gaussian to a quasi-Lorentzian 
line shape may be driven by an increasing dilution of 
(nearly) static magnetic impurities or defects, and vice 
versa. Thus, the larger disorder probed by T 2 ~ 4 peaks 
consistently coincides with a larger concentration of mag- 
netic defects indicated by (3 > 1, which agrees with the 
shortening of the magnetic coherence length observed on 
approaching the field-induced transitions^ 

In contrast, spin-lattice relaxations are governed by in- 
dependent excitations, whose effect on T2 is unimportant 
as it is masked by the relaxation channel sketched above. 
The thermally activated behavior of reveals a gap 
A of order 100 K in their spectrum. A thorough investi- 
gation of A as a function of B ext for the different chain 
types, and a fitting of the experimental energy gaps to 
a theoretical model, are the subject of a separate paper. 
Here we anticipate that, according to the model, the ac- 
tivated spin-lattice relaxations are essentially driven by 
the spin-flip of a single Co 3+ moment. 

D. Spin freezing: relevance for the magnetization 
steps 

We conclude this discussion with a comment on the 
low-temperature multi-step behavior in M(H). The 
magnetic structures underlying each step could not be ac- 
cessed in our NMR study, due to the non-spherical shape 
of our sample, inducing demagnetization-dependent line 
broadening and offsets which prevent finer resolution of 
dipolar fields, as well as to exceedingly slow spin-lattice 
relaxations below 5 K. Nevertheless, T2 data at 7-15 K re- 
veal a scenario of glassy spin dynamics involving freezing 
of sizable magnetic aggregates and large quasi-static in- 
homogeneities close to the field-induced transitions. The 
collective arrangement of such frozen entities at lower 
temperature is unknown. We may argue however that 
it gives rise to complex field-dependent micromagnetic 
configurations or superstructures, resulting in a multi- 
plicity of local minima separated by potential barriers in 
the free energy, becoming metastable below 2 K. If this 
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picture holds true, however, the physical mechanism be- 
hind magnetization steps should be understood in terms 
of thermally assisted hopping between essentially meso- 
scopic metastable configurations, probably incommensu- 
rate with the crystal lattice, rather than quantum tun- 
neling as in molecular magnets. 
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Appendix A: The ° Co resonance 

The resonance frequencies of 59 Co at the non-magnetic 
site I, for an arbitrary experimental geometry, depend 
on the external field and sample alignment in a rather 
complicated way, which deserves some comment. 

The field at the nucleus B nuc (proportional to the 
resonance frequency vz = 7-B n „ c ) differs from the lo- 
cal field Bi oc at the ion site by the chemical shift, 
B nuc = (1 + K)Bi oc , where the chemical shift tensor K 
is usually large and anisotropic in cobalt^ In a magnetic 
material, the local field is the vector composition of the 



external and the internal field, Bi, 



B ext +B 



where 



the internal field Bi nt originates from the electronic mo- 
ments. The overall dependence of vz is summarized by 
Eqn.El 

In the present case, the sources of -Bj„t are i) the dipo- 
lar coupling with neighboring Coll spins; ii) the net de- 
magnetization field, namely, the difference between the 
demagnetization effects of the outer sample surface and 
a Lorcntz sphere, which does not vanish due to the non- 
spherical shape of the sample; and Hi) possibly, a trans- 
ferred hyperfine contribution, due to the polarization of 
wave functions. Even in a hard easy-axis magnet like 
the present one, the direction of magnetic moments with 
respect to the crystal axes, hence of Bi nt , depends on 
H fa Hq Bext via the finite transverse magnetic sus- 
ceptibility Xab of the system, leading to a field-induced 
tilting of the Co 3+ spins by a small but non-negligible an- 
gle ip ex XabH±, and by w —'0/2 for Bi nt , if the dipolar 
contribution dominates. 

In a crystal, 59 Co (I = 7/2) is also coupled by its 
quadrupole moment to the local electric field gradient 
(EFG) via a quadrupolar spin Hamiltonian of the form 



H 



Q 



6 



Ml -1(1+1) 



n 



H + r 



(Al) 



where I± = I x ± il y , and I x I y I z are the nuclear spin 
components in the crystal reference frame. Here the 
quadrupolar frequency vq is defined as a function of the 
nuclear quadrupole moment Q and the main component 
V zz of the EFG tensor as v Q = 3eV zz Q/[2hI(2I-l)], and 



In the presence of a large magnetic field (either ex- 
ternal or internal) at the nucleus, such that Vz ~3> vq, 
the quadrupolar interaction of Eqn. lAll bchaves as a per- 
turbation, and splits the nuclear Zeeman transitions into 
a multiplct of 21 satellite lines. In a second-order ex- 
pansion vs. vq/vz, the frequency v n of the \n— 1/2) 
|n+l/2) Zeeman transition (n = -27 + 1/2,-27 + 
3/2, ... 2/ — 1/2) depends on B nuc as 



H 

vz 



(7+1/2) 2 -! [2S(0,0)-C(0,0)] 



v Q A(9, <j>)n 



vz 



[4B(0,(t))-C{9,(f))}n 2 (A2) 



where 9, </> are the polar coordinates of B nuc in the 
crystal frame, and the coefficients A, B, C are defined as 

A(9,(/)) = i (3 cos 2 (9- 1 + 77sin 2 6»cos2(/)) 

B{9, <P) = — sin 2 9 [cos 2 9 (9 - 6i] cos 2</> + rj 2 cos 2 2<j>) 
+ if sm 2 2(f>) 

C(9,cj)) = -^{9sin 4 6> + 6?7sm 2 6>(l + cos 2 6l)cos2(/> + 

rj 2 [cos 2 20 + 2 cos 2 9 (l + sin 2 2</>) + cos 2 9 cos 2 20] } 

(A3) 

We note that all second order terms vanish for a cylin- 
drical EFG (rj = 0) and EFG principal axis collinear to 
Bnuc (0 = 0). 



Appendix B: The spin-lattice relaxation function 

The approach to thermodynamic equilibrium of an 
ensemble of nuclear spins / is governed by a closed 
set of rate equations for the m-th level population N m 
(m = —I)- In the presence of a magnetic relax- 

ation channel, it has the formic 



N„ 



\Vyy \ ~ \V X x\\/V zz is the EFG rhombicity factor 



34 



= W[(I + m + l)(I-m)(N m+1 -N^ +1 ) 
+ (I - m + l)(I + m)(N m ^ ^ N^l,) 
- 2(I 2 - m 2 + I)(N m -N^)] (Bl) 



where is the population at equilibrium, and W is 

the transition probability \m) -f-> \m — 1) between two 
adjacent Zeeman levels. In the case of quadrupolc-split 
NMR lines with the selective excitation of a single nu- 
clear transition of order n (n = I — 1/2, . . . , —I + 1/2), 
the signal intensity is proportional to the population dif- 
ference A n = N n+1 / 2 —N n+1 /2 of the two levels involved, 
rather than the expectation value of the total spin (I z ). 
Solution of Eqn. IB 1 1 yields a recovery of A n according 
to the superposition of 21 exponential components for 
generic n (reduced to 7+1/2 for n = 0) with relative 
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weights depending on the saturation method as well as 
on n. 

In our experiments, the spin ensemble was prepared 
with a fast saturation, i.e. with a pulse train much shorter 
than the spin-lattice relaxation time T\ = (2W)^ 1 , so 
that populations N m of nuclear levels m ^ n ± 1/2 not 
involved in the transition were practically unaffected . 42 ' 43 
For the central line n = 0ofa/ = 7/2 nuclear species as 
59 Co, the recovery law following such an initial condition 



is then calculated as follows. 
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